% Generate distribution of firm efficiency from standard normal draws U and parameters

function [phi] = DrawProductivity(U,param)
lnphi = param.std_phi*U; % log efficiency given desired standard deviation
phi = exp(lnphi);
phi = phi / mean(phi); % Ensure that efficiency has mean one
end
